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Abstract 

We study the stability of a scalar inflaton field and analyze its point 
attractors in the phase space. We show that the value of the inflaton field in 
the vacuum is a bifurcation parameter and prove the possible existence of a 
limit cycle by using analytical and numerical arguments. 
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The phenomenon of chaos^, which has become very popular, rejuvenated 
interest in nonUnear dynamics, and the chaotic behaviour of classical field 
theories is currently subject of intensive research^'^'^'^. In this respect it is 
of great interest to investigate the existence of some precursor phenomena 
of chaotic motion. In particular nonlinear effects should lead to bifurcations 
which qualitatively change the system properties^'^'^. 

In this paper we study the stabihty of a scalar inflaton field^'^°'^^ and 
analyze its bifurcation properties in the framework of the dynamical system 
theory^'''''^. The model we study is very schematic, so it can be seen as a 
toy model for classical nonlinear dynamics, with the attractive feature that 
it emerges form particle physics and cosmology. 

It is generally believed that the universe, at a very early stage after the big 
bang, exhibited a short period of exponential expansion, the so-called infla- 
tionary phase^'^°'^^. In fact the assumption of an inflationary universe solves 
three major cosmological problems: the flatness problem, the homogeneity 
problem, and the formation of structure problem. 

The Priedmann-Robertson- Walker metric^^ of a homogeneous and isotropic 
expanding universe is given by: 

ds" = de - a^it)\ , , + r'^ide'' + sin^ edip")], (1) 
1 — kr^ 

where k — 1, —1, or for a closed, open, or flat universe, and a{t) is the scale 
factor of the universe. 

The evolution of the scale factor a{t) is given by the Einstein equations: 

d = — —G{p + 3p)a, 
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where p is the energy density of matter in the universe, and p its pressure. 
The gravitational constant G = (with h = c = 1), where Mp = 1.2 ■ lO"*^^ 
GeV is the Plank mass, and H — a/ a is the Hubble "constant", which in 
general is a function of time. 

All inflationary models postulate the existence of a scalar field 0, the 
so-called inflation field, with Lagrangian^^: 

L = ^d^<i>d^<i> - y (0) (3) 

where the potential V{(j)) depends on the type of inflation model considered. 
Here we choose a real field but also complex scalar can be used^^ . The scalar 
field, if minimally coupled to gravity, satisfies the equation^'^°'^^: 

□0 = + 3(-)0-lv20 = -^, (4) 

where □ is the covariant d'Alembertian operator. We suppose that in the 
universe there is only the infiaton field, so the Hubble "constant" H is related 
to the energy density of the field by: 



,,^4 = (V4 = ?^[f + (M! + vm (5) 
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Immediately after the onset of inflation, the cosmological scale factor 
grows exponentially^'^^'^-*^. Thus the term V^0/a^ is generally believed to be 
negligible and, if the infiaton field is sufficiently uniform (i.e. 0^, (V0)^ << 
y(0)), we end up with a classical nonlinear scalar field theory in one dimen- 
sion: 

dV 

+ m{<^)<^ + ^ = 0, (6) 



where the Hubble " constant" H is an exphcit function of 0: 

= ^VW. (7) 

Most inflation scenarios are just formulated in this classical setting. In fact 
the quantization of the inflation scenario can be regarded to be still an open 
problem^^. 

The classical infiaton dynamics can be studied in the framework of the 
dynamical system theory^'^'^. The second order equation of motion can be 
written as a system of two first order differential equations: 

x^g{x) (8) 

where x = (0, x) is a point in the two dimensional phase space and g — 
(91:92) is given by: 

^i(0,x)=X, ^2(0,x) = -3//(0)x-^^. (9) 
The system is non-conservative because the function 

d^v{9) = ^ + ^ = -m<P) (10) 

is not identically zero. The fixed points of the system are those for which 
9i{(p, X) = and t/2(0, x) = 0, i.e: 

The deviation 6x{t) = x{t) — x{t) from the two initially neighboring tra- 
jectories X and X in the phase space satisfies the linearized equations of 
motion: 

^6x{t) = T{t)Sx (12) 



where r{t) is the stabihty matrix: 

■^«^(_^°3,|| -si,) I- 

At least if an eigenvalue of r(t) is real the separation of the trajectories grows 
exponentially and the motion is unstable. Imaginary eigenvalues correspond 
to stable motion. In the limit of time that goes to infinity the eigenvalues 
of the stability matrix are the Lyapunov exponent^. It is well known that 
for two-dimensional dynamical system the Lyapunov exponent can not be 
positive and so the system is not chaotic, i.e. there is not global instability. 

However, we can be assured that the universe is crowded with many 
interacting fields of which the inflaton is but one^^. The nonlinear nature of 
these interactions can result in a complex chaotic evolution of the universe 
and the local instability of the inflaton field is a precursor phenomenon of 
chaotic motion. Following the Toda criterion^^, we assume that the time 
dependence can be eliminated, i.e. r(0(t),x(t)) = r(0, x). 

The eigenvalues of the stability matrix are given by: 



S 1 / r)2T/ f)fr 

<'M = -2i^W±2^9i^n^)-4g^-12Xg^. (14) 

The pair of eigenvalues become real and there is exponential separation of 

neighboring trajectories, i.e. unstable motion, if: 

d'^V dH , , 

+ 3x^<0. (15) 



Particularly when x = 0, e.g. the fixed points, we obtain local instability 
when: 

IS < 0. (la) 



i.e. for negative curvature of the potential energy. The fixed points are stable 
if they are point of local minimum of V^(0) and unstable if are points of local 
maximum. 

As stressed at the beginning the potential depends on the type 

of inflation model considered, and it is usually some kind of double-well 
potential. We choose a symmetric double-well potential: 

= (17) 

where ±v are the values of the inflaton field in the vacuum, i.e. the points 
of minimal energy of the system. 

We observe that the inflaton field value in the vacuum v is a bifurcation 
parameter. Bifurcation is used to indicate a qualitative change in the features 
of the system under the variation of one or more parameters on which the 
system depends^'^'^. First of all we consider the case v = 0, i.e. V{(j)) = 
(A/4)0^. In this situation there is only one fixed point {(f)* — 0, x* — 0) 
which is a stable one being: 

9^ - 3A0^ > 0. (18) 

The fixed point (0* = 0, x* = 0) is a point attractor, as shown in Fig. 1 
where we use a fourth order Runge-Kutta method^^ to compute the classical 
trajectories in the phase space. 

Instead ior v ^ there are three fixed points: 

(0* = O,x* = O), (0*=^,x* = O), (0* = -^;,x* = O), (19) 

and the condition for the instability becomes: 

V V 



Obviously (0* = 0, %* = 0) is an unstable fixed point, and in particular a 
saddle point because the stability matrix has real and opposite eigenvalues. 
On the other hand (0* = ±i>, x* — 0) are stable fixed points. 

Form eq. (7) we find four possible functions for the Hubble "constant": 

H{4>)^±^\4>'-v\ (21) 

but also: 

i^(0) = ±7(0'-^'), (22) 

where 7 — ^J2^^GX/3 is the dissipation parameter. The choice of the Hubble 
function is crucial for the dynamical evolution of the system. 

In certain non-conservative systems we could find closed trajectories or 
limit cycles toward which the neighboring trajectories spiral on both sides. 
It is sometime possible to know that no limit cycle exist and the Bendixson 
criterion^^, which establishes a condition for the non-existence of closed tra- 
jectories, is useful in some cases. Bendixson criterion is as follows: if div{g) 
is not zero and does not change its sign within a domain D of the phase 
space, no closed trajectories can exist in that domain. 

In our case we have div{g) — — 3if(0), and so the presence of periodic 
orbit is related to the sign of H{(p). If H{(f>) = 7|0^ — u^l we do not find 
periodic orbits and the inflaton field goes to one of its two stable fixed points, 
which are points attractors (see Fig. 2). The vacuum is degenerate but if 
we choose an initial condition around the saddle point there is a dynamical 
symmetry breaking^^ towards the positive v or negative —v value of the 
infiaton field in the vacuum. This symmetry breaking is unstable because 
neighbour initial conditions can go in different point attractors. 
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Instead, if we choose H{(j)) = 7(0^ — v^) the numerical calculations of 
Fig. 3 show that exists a limit cycle, the two stable fixed points are not 
point attractors, and the inflaton field oscillates forever. Obviously more 
large is v more large is the limit cycle. 

In conclusion, we have studied the stability of a scalar inflaton field. 
With a symmetric double-well self energy the value of the inflaton field in 
the vacuum is a bifurcation parameter which changes dramatically the phase 
space structure. The main point is that for some functional solutions of the 
Hubble "constant" the system goes to a hmit cycle, i.e. to a periodic orbit. 

The inflaton field is not chaotic but its local instability can give rise to 
a complex chaotic evolution of the universe due to its nonlinear interactions 
with other fields. In the future it will be very interesting to study these 
effects which can perhaps lead to some observable implications like a fractal 
pattern in the spectrum of density fluctuations. 

The author is greatly indebted to Prof. S. Graffi, Prof. V. R. Manfredi 
and Dr. A. Riotto for many enlightening discussions. The author acknowl- 
edges Prof. J. M. G. Gomez for his hospitality at the Department of Atomic, 
Molecular and Nuclear Physics of " Complutense" University, and the "Ing. 
Aldo Gini" Foundation for a partial support. 
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Figure Captions 



Figure 1: Inflaton field time (up), and its phase space trajectory (down), 
with 7 = 1/2, A = 1 and v = 0. Initial conditions: = and X = = 1/2. 
Figure 2: Inflaton field vs time (up), and its phase space trajectory (down), 
for -ff(0) = 7|0^ — v^l with 7 = 1/2, A = 1 and v — 1. Initial conditions: 
(/) = 0andx = </' = l/2. 

Figure 3: Inflaton field vs time (up), and its phase space trajectory (down), 
for i/(0) = 7(0^ — v^) with 7 = 1/2, A = 1 and v = 1. Initial conditions: 
= and x = = 1/2. 
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